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ABSTRACT 

The potential of tidal disruption of stars to probe otherwise quiescent supermassive 
black holes cannot be exploited, if their dynamics is not fully understood. So far, 
the observational appearance of these events has been derived from analytical ex¬ 
trapolations of the debris dynamical properties just after disruption. By means of 
hydrodynamical simulations, we investigate the subsequent fallback of the stream of 
debris towards the black hole for stars already bound to the black hole on eccentric 
orbits. We demonstrate that the debris circularize due to relativistic apsidal precession 
which causes the stream to self-cross. The circularization timescale varies between 1 
and 10 times the period of the star, being shorter for more eccentric and/or deeper 
encounters. This self-crossing leads to the formation of shocks that increase the ther¬ 
mal energy of the debris. If this thermal energy is efficiently radiated away, the debris 
settle in a narrow ring at the circularization radius with shock-induced luminosities 
of - 10 - 10 3 LecM- If instead cooling is impeded, the debris form an extended torus 
located between the circularization radius and the semi-major axis of the star with 
heating rates ^ 1 — 10 2 I/Edd- Extrapolating our results to parabolic orbits, we infer 
that circularization would occur via the same mechanism in ^ 1 period of the most 
bound debris for deeply penetrating encounters to ^ 10 for grazing ones. We also 
anticipate the same effect of the cooling efficiency on the structure of the disc with 
associated luminosities of ~ 1 — lO^Edd an( ! heating rates of ~ 0.1 — ll^Edd- In the 
latter case of inefficient cooling, we deduce a viscous timescale generally shorter than 
the circularization timescale. This suggests an accretion rate through the disc tracing 
the fallback rate, if viscosity starts acting promptly. 

Key words: black hole physics - accretion discs - hydrodynamics - galaxies: nuclei. 


1 INTRODUCTION 

Most galaxies have been found to contain a supermassive 
black hole (SMBH) at their centre orbited by stars. If one 
of these stars wanders within the tidal radius of the black 
hole, the tidal force of the black hole exceeds the star’s self¬ 
gravity and the star is torn apart. Such an event is called 
a tidal disruption event (TDE). After the disruption, the 
stellar debris evolve to form an elongated stream of gas that 
falls back towards the black hole. In the standard picture of 
TDEs, these debris then circularize to form an accretion disc 
where they accrete viscously e mitting a thermal flare mainly 
in the UV to soft X-ray band (|Lodato fc Rossill201lh . This 
flare could also be accompanied by a radio signal associ¬ 
ated to a relativistic jet originating from the inner region 


* E-mail: bonnerot@strw.leidenuniv.nl 


of the disc (|Giannios &; Metzgerl 1201 li b A handful of candi¬ 
date TDEs have been detected so far in these bands (e.g . 
Esauei et al. 200p: iKomossa et al 1 120041: iGezari et al.l 12009 : 


van Velzen et al.ll201ll: IZauderer et al.l 1201 lh~ 


TDEs are powerful tools to detect SMBHs in other¬ 
wise quiescent galaxies. Furthermore, they can in principle 
be used to estimate the properties of both the black hole 
and the disrupted star and to probe the physics of accre¬ 
tion and relativistic jets. In practice, deriving such con¬ 
straints from observations is challenging as it requires a 
precise understanding of the dynamics of TDEs. The lat¬ 
ter has been the focus of many analytical and numerical in¬ 
vestigations undertaken sin ce the 80s. A dist i nctive featur e 
of the pioneering w o rks bylLacv et al.l dl982h. Reesl (Jl988), 
lEvans fc Kochanekl (Il989h and Phinnevl (|l989h is a t _5//3 
decrease of the rate at which the stellar debris fall back 
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towards the black hole0 The same sloi 
ted to the observed TDE light curves ( 

pe is generally fit- 

Esauei et al.l 20081: 

Gezari et al.ll2008l. 20091: ICaDoelluti et a 

. 2009lj although it 

is only expected in the X-ray band (Loc 

ato & Rossi! 1201 lh . 

It assumes that the accretion rate, and therefore the lumi- 


nosity, traces the fallback rate of the debris. In turn, this 
requires that the debris circularize to form a disc and that 
this disc is accreted faster than it is fed from the fallback 
stream. However, the mechanism driving the circularization 
is still unknown. 

Various effects are likely to be involved in this pro¬ 
cess, whose common feature is to dissipate the kinetic en¬ 
ergy of the debris, injecting a large amount of thermal en¬ 
ergy into the newly formed disc. However, the efficiency of 
this energy transfer is certainly dependent on the mech- 
anism considered. The main candidates are the following 
([Evans V Kochaneklll989l : iKochanekl 1 1 994| j : 


(i) Pancake shock : As the star is disrupted, part of its 
material is accelerated out of the initial orbital plane. As 
a result, the debris inside the stream have a range of incli¬ 
nations. Their orbits are therefore vertically focussed and 
intersect the orbital plane close to pericentre. At this point, 
the stream is strongly compressed, leading to the formation 
of a pancake shock. 

(ii) Self-crossing: When they reach pericentre, the debris 
experience changes of their apsidal angles driven by rela¬ 
tivistic apsidal precession or hydrodynamical effects. This 
can lead to the self-crossing of the stream: as the leading 
parts move away from the black hole after pericentre pas¬ 
sage, they collide with the part that is still falling back, 
generating shocks. 

(iii) Shearing : As the stream comes back to pericentre, 
the orbits of the debris are radially focussed due to their 
large range of apocentres but small range of pericentres. 
This effect is similar to a passage into a nozzle. The debris 
experience shearing at this point as they have a range of 
apsidal angles and eccentricities induced by the disruption 
of the star. This effect is enhanced by relativistic apsidal 
precession: as they have a range of pericentre distances, the 
orbits of the debris precess by different angles, leading to 
further shearing. 


If an accretion disc forms from the debris, its structure 
and evolution are two additional uncertainties in the mod¬ 
els of TD Es. They depend on the re lative efficiency of three 
processes (lEvans fe Kochaneklll989l j: circularization, viscous 
accretion and radiative cooling. Denoting the timescales of 
these processes by Girc, t w i SC and t C ooi respectively, three lim¬ 
iting regimes are to be expected. In the case t v i SC < t c i rc , 
viscosity may be important during the circularization pro¬ 
cess. If instead t v i SC > t c irc, accretion begins only once a disc 
is formed. In this case, if t coo \ < t c i rc , the disc cools during 
its formation and is therefore geometrically thin. If instead 
tdrc < tcooi, the disc puffs up while it forms due to its ex¬ 
cess of thermal energy. Many author s have assumed that 
the d is c is geometrically th i n or sl i m (IStrubbe V Quataert 
20091: IC annizzo V Gehr elsl 200 91: Cannizz o et al.l l201ll ; 


Lodato V Rossill201ll : IShen Matznerll2014h using~the stan 


1 Th is rate was later shown to be dependent on the stellar struc- 
ture (iLodato et alJl2009t iGuillochon fe Ramirez-Ruidl2013h . 


dard a parametrization (IShakura V Sunvaevlll973 j. Other 
investigations considered the poss i bility of a geometri¬ 
cally thick disc (iLoeb V Ulmer|[l997l : I Coughlin V Begelmanl 
12014 . 

Numerical simulati o ns of tidal_disruptions 

(lEvans V Kochanekl Il989l : ILodato et al.l l2009h have of¬ 
ten used smoothed particle hydrodynamics (SPH) primarily 
because of its ability to deal with large regions of space 
devoid of gas. This technique is also well suited to simulate 
the subsequent fallback of the debris towards the black 
hole. The computational cost of such a simulation scales 
with the total number of SPH particles used to model the 
stream. To accurately follow their evolution, each part of 
the stream must contain enough particles. This condition is 
hard to fulfill for a long stream because the SPH particles 
are spread on a large volume. The length of the stream 
increases with the stellar to black hole mass ratio q and 
with the eccentricity e of the star. The typical values of 
these parameters are q = 10 6 and e — 1. This extreme 
value of e comes from the fact that most disrupted stars 
are scattered from the sphere of influence of the SMBH, 
which is much larger than the tidal radius (iFrank fe Reesl 
Il976h . For the Milky Way, the ratio of th e se tw o distances 
is ~ 10 5 . However, as noted bv lAval et al.l (200 0), following 
the fallback of the debris for these typical values of q and e 
numerically is a computational burden. In their simulation, 
the leading part of the stream is composed of very few 
SPH particles that come back almost one by one towards 
the black hole. As a result, the evolution of these particles 
cannot be followed accurately. Even if they use a relatively 
low number of particles (N = 4295), this issue is so extreme 
that it is likely to persist for larger particle numbers. The 
high computational cost of such a simulation is not limited 
to SPH but generalizes to other computational techniques. 

In the few other investigations of this problem, either q 
or e has therefore been lo wered. The first option has been 
chosen by several auth ors (Ra mirez-Ruiz fe Rosswo3 f2009: 
IGuillochon et al.l l2014j who investigated q = 10 3 and e = 1. 
The physical motivation in this case is the tidal disruption 
of a star by an intermediate mass black hole. Both pancake 
shock and shearing were shown to be inefficient at circular¬ 
izing the debris in this case. Instead, the main effect seen in 
these simulations is the expansion of the stream caused by 
its pericentre passage. Circularization is more likely driven 
by self-crossing either indu ced b y this expansion or by rela¬ 
tivistic apsidal precession. lHavasaki et al.l (|2013 ;j made the 
second choice, as in the present paper, which corresponds to 
tidal disruptions of bound stars. They considered q = 10 6 
and e = 0.8 and found that self-crossing driven by rela¬ 
tivistic precession efficiently drives the circularization of the 
debris in this case. 

Several mechanisms are able to put a star on a bound 
orbit entering the tidal radius of a SMBH. The first involves 
a unequal mass SMBH binary interacting with a surrounding 
stellar cusp. In this situation, bound stars can be scattered 
inside the tidal radius of the primary black hole through 
a combination of Kozai inter actions and close encou nters 
with the secondary black hole ([Chen et al.ll2009l . 1201 lh . The 
second mechanism takes place during the coalescence of a 
SMBH binary. Due to the anisotropic emission of gravita¬ 
tional waves in this phase, the remnant black hole undergoes 
a kick. It can lead to close encounters between the black 
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hole and the surrounding sta rs leading to their tidal dis¬ 
ruptions (IStone &; Loebl l201lh . Furthermore, mean-motion 
resonances can pull surrounding s tars inwards during the 
final phase of the SMBH merger (ISeto &; Mutol 2010) en¬ 
hancing the rate of tidal disruptions of bound stars by this 
mechanism. The third mechanism takes place after the tidal 
separation of a stellar binary by a black hole, which places 
one of the components on an eccentric orbit. Due to both 
encounters with the other stars orbiting the SMBH and grav¬ 
itational wave emission, its orbit shrinks and circularize. Un¬ 
der certain conditions, it can t hen be tidally disru pted by 
the SMBH on a bound orbit (|Amaro-Seoane et al.ll2Q12h . 

In this paper, we use SPH simulations to investigate the 
circularization process leading to the formation of the disc. 
We also characterize the structure and evolution of this disc. 
Relativistic apsidal precession, involved in the circulariza¬ 
tion proces^j^treate^exa^tlyv^iUjnosWprevmus stud¬ 
ies (I Ram irez-Ru iz &; Rosswod 120091 : iGuillochon et al H'2014 
iHavasakT^Tlir 1201.'4 used an approximate treatment of this 
effect. The parameters that are likely to play a significant 
role in the circularizaton process have also been varied. Two 
extreme cooling efficiencies, encoded in the equation of state 
(EOS) of the gas, have been considered. The effect of the or¬ 
bit of the star has also been examined, by varying both its 
penetration factor and its eccentricity. Our results demon¬ 
strate that circularization is a fast process, happening in a 
few orbits of the debris around the black hole, and driven 
mostly by relativistic apsidal precession. They also confirm 
the expected effect of the cooling efficiency on the struc¬ 
ture of the resulting disc. In addition, we found different 
channels of circularization for different orbits of the star. In 
particular, by increasing its eccentricity, we got insight into 
the circularization process at work in the standard parabolic 
case. 

This paper is organized as follows. In Section [21 we de¬ 
scribe the SPH simulations performed. The results of these 
simulations are presented in Section [3] varying the follow¬ 
ing parameters: the gravitational potential, the EOS, the 
depth of the encounter and the eccentricity of the star. The 
effect of the resolution on these results is also examined. 
These results ar e discussed in Section 5] a nd compared to 
other studies by IShiokawa et al.l (|2015h and lHavasaki et al.l 
(l2015l l that appeared at the same time as this paper. 


2 SPH SIMULATIONS 

We investigate the tidal disruption of a star of mass M* = 
1 M 0 and radius R* = 1 R© by a non-rotating black hole 
of mass Mh = 10 6 M©. In this configuration, the tidal ra¬ 
dius is Rt = R*(M\ 1 /M*) 1 / 3 = 100 R©. Several initial ellip¬ 
tical orbits of the star are considered. They have different 
pericentre distances R p defined via the penetration factor 
/ 3 = Rt/Rp by setting P — 1 and /3 = 5. It corresponds to 
pericentre distances R p = 100 R© and R p = 20 R© respec¬ 
tively. For p = 1, only an eccentricity e = 0.8 is considered. 
For p — 5, two different values e = 0.8 and e = 0.95 are 
investigated. For these orbits, all the debris produced by 
the disruption stay bound to the black hole. This is gener¬ 
ally the case if e < e cr it = 1 — (2//3)(Mh/M*) -1 / 3 , where 
e C rit = 0.996 and 0.98 for p = 5 and p = 1 respectively. The 
semi-major axis a* = R P /( 1 — e) of these orbits ranges be- 


- SPH simulation 

-- Relativistic potential 

. Schwarzschild metric 



R t 


Figure 1 . Trajectory of the centre of mass of the gas in the simu¬ 
lation corresponding to model RI5e.8 (black solid line) compared 
to that of a test particle on the same orbit in the relativistic 
potential (red dashed line) and in the Schwarzschild metric (blue 
dotted line). All the trajectories start at the initial position of the 
star indicated by the brown point on the right of the figure. The 
black arrow specifies the direction of motion. The trajectory of the 
test particle in the relativistic potential and in the Schwarzschild 
metric is followed for 4 orbits. That of the centre of mass of the 
gas in the simulation is followed during the disruption of the star 
and the fallback of the debris until t/P* = 8. The transition be¬ 
tween the disruption and the fallback phase, treated with two 
different codes, is indicated by the green dash perpendicular to 
this trajectory. After a few pericentre passages, the motion of the 
gas is affected by hydrodynamical effects and the trajectory of its 
centre of mass differs from that of a test particle. 


tween 100 R© and 500 R©, corresponding to orbital periods 
P* = 27 t (GMh/a*) 1 ^ 2 between 2.8 h and 31 h. 

Both the disruption of the star and the subsequent fall¬ 
back of the debris towards the disruption site are simulated 
using SPH. To increase efficiency, the simulations of these 
two phases are performed separately and make use of two 
different codes, each adapted to a specific phase. For the 
disruption phase, we use a code that takes into account 
self-gravity ([Bate et al.lll995h . For the fallback phase, we do 
not consider self-g ravity and make use of the highly efficient 
code PHANTOM (jPrice & Federrat3l2Qld : lLodato & Pricel 
l201ftl that is optimized for studying non-self-gravitating 
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problem^. This choice is legitimate because the self-gravity 
force is only needed during the disruption phase where it op¬ 
poses the tidal force of the black hole. In the fallback phase, 
the gravitational interactions between the debris is negligi¬ 
ble compared to their hydrodynamical interactions and the 
tidal force of the black hole. The disruption phase is followed 
until the most bound debris come back to pericentre. Their 
properties are recorded at this point and constitute the ini¬ 
tial conditions for the simulation of the fallback phase. As 
this paper’s primary focus is on the fallback of the debris, 
the disruption phase is only simulated to get these initial 
conditions and not discussed in detail. 

The co de used for th e disru ption phase has already been 
adopted bv lLodato et al.l (120091 ) to simulate tidal disruptions 
of stars on parabolic orbits. In the present paper, we follow 
the same procedure and numerical setup but consider ellip¬ 
tical orbits instead. Also using the same method, we model 
the star as a polytropic sphere with 7 = 5/3 containing 
100K particles. However, as explained in subsection 13.41 the 
fallback phase is simulated at a higher resolution with the 
stream of debris containing 500K particles. 

In order to resolve the shocks, the code used for the 
fallback phase, PHANTOM, includes the standard artifi¬ 
cial viscosity prescription that depends on the parameters 
a AV and /3 AV . In addition, the lMorris & Monaghanl (Il997l ) 
switch is implemented to reduce artificial dissipation away 
from shocks. To this end, a AV is allowed to vary between 
two values and according to a source and decay 

equation. In this paper, we use a A ^ = 0.01, a A ^ x = 1 and 

p AV = 2 . 

The orbits considered for the star have pericentre dis¬ 
tances comparable with the gravitational radius R g = 
GM^/c 2 of the black hole. More precisely, R p = 47 R g and 
R p — 9AR g for /3 — 1 and [3 — 5. Therefore, relativistic 
effects must significantly affect the motion of the gas when 
it passes at pericentre. One of these effects is relativistic ap- 
sidal precession, a mechanism involved in the circularization 
process causing the orbit of a test particle to precess at each 
pericentre passage by a given precession angle. For the or¬ 
bits considered, the precession angle varies between 13.5 and 
89.7 degrees for p = 1 and /3 = 5 respectively. In order to 
investigate the influence of relativistic effects, the black hole 
is modelled by an external potential that is either Keplerian 
or relativistic both in the disruption and the fallback phase. 
In the following, the variables are labelled by the letter K 
if they are computed in the Keplerian potential and by the 
letter R if they are computed in the relativistic potential. 
These potentials are respectively given by 


<f> K 

<f> R 


GM h 
R ’ 

GAP _ ( 2 \ r f R-Rz \ 2 vj_ 

R \R-2R S ) [\R-2R g J 1 2 


( 1 ) 

( 2 ) 


where v T and v t are the radial and tangential velocity of 
a test particle respe ctively. We adopt this re lativistic po¬ 
tential, designed by iTeieda fc Rosswoel ([2013 ), because it 
contains an exact treatment of apsidal precession around 
a non-rotating black hole. In other words, the value of 


2 Self-gravity is now implemented into PHANTOM. However, it 
was not available at the beginning of this work. 


the precession angle for a test particle on a given orbit is 
the same in this potential and in the Schwarzschild met¬ 
ric. This is an improvement compared to most previous in¬ 
vest i gations that used an approximate treatment of this ef¬ 
fect (iRa mirez-Ruiz fc Rosswodl2009l : iGuillochon et al.1l2014l : 
lHavasaki et al.l 20131 ). 

For the disruption phase, the star is initially placed at 
the apocentre of its orbit, a distance R a = (1 + e)R p /(l — e) 
from the black hole. Its initial velocity v a is determined by 
conservation of specific orbital energy and angular momen¬ 
tum between R a and R p . In the Keplerian and relativistic 
potentials, the specific orbital energy is given by 


1 

2 


(Vr +Vt) 


GM h 

R 


1 / R 

2 \R-2R g 


2 


2 




+ 


R \ 2 

R — 2R g ) Vt 


(3) 


GMh 

R 


(4) 


and the specific angular momentum by 

l K = Rv t , 

R _ Rv t 

R-2R g 


(5) 

(6) 


respectively. The initial velocity of the star at apocentre is 
therefore 



R = / GMh \ 1/2 2^R 2 {R 1 -2R g ) _ 

a V .Ra J RR(R a + Rp)(Rp-‘2R e )+2R s Rl) 1 / 2 ’ ( ' 


depending on the potential. 

In order to demonstrate the correct implementation of 
the relativistic potential into the SPH codes, we anticipate 
the next section and analyse the motion of the gas in the 
simulation associated to model RI5e.8. In this model, the 
relativistic potential is used and the orbit of the star has a 
penetration factor = 5 and an eccentricity e = 0.8. Fig. 
\T\ shows the trajectory of the centre of mass of the gas in 
the simulation corresponding to model RI5e.8 compared to 
that of a test particle on the same orbit in the relativistic 
potential and in the Schwarzschild metric. The trajectory 
of the test particle in the relativistic potential and in the 
Schwarzschild metric is followed for 4 orbits. That of the 
centre of mass of the gas is followed during both the disrup¬ 
tion of the star and the fallback of the debris until t/P* = 8 . 
The trajectory of the test particle is the same in the rela¬ 
tivistic potential and in the Schwarzchild metric, confirming 
that apsidal precession is treated accurately i n the relativis- 
tic potential as found bv ITeieda fe Rosswod ([20131 ). During 
the disruption and the beginning of the fallback phase, the 
centre of mass of the gas and the test particle have the same 
trajectory, validating the implementation of the potential 
into the SPH codes. After a few pericentre passages, the 
two trajectories start to differ as the hydrodynamical effects 
on the debris become dominant. This evolution will be in¬ 
vestigated in detail in the next section. 

The accretion onto the black hole is modelled by an 
accretion radius fixed at the innermost stable circular orbit 
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Figure 2. Snapshots of the fallback of the debris at different times £/P* = 0, 1.2,3 and 8 for models RI5e.8 (upper panel) and KI5e.8 
(lower panel). For these models, the period of the star P* is 2.8 h. The colours correspond to the column density E of the gas whose 
value is indicated on the colour bar. The white point represents the black hole. The dashed white circle on the last snapshot represents 
the circularization radius given by equations ca and m for the Keplerian and relativistic potentials respectively. 


of 6P g for a non-rotating black hole. Particles entering this 
radius are removed from the simulations. 


The simulation of the fallback phase is performed for 
two different equations of state (EOS) for the gas: locally 
isothermal and adiabatic. For the locally isothermal EOS, 
each SPH particle keeps their initial specific thermal ener¬ 
gies. Physically, these two EOS represent the two extreme 
cases for the rate at which an excess of thermal energy is 
radiated away from the gas. For the locally isothermal one, 
every increase of thermal energy is instantaneously radiated 
away while for the adiabatic one, none of this energy is ra¬ 
diated. 


Four parameters are therefore considered to simulate 
the fallback of the debris: the potential, the EOS, the pen¬ 
etration factor p and the eccentricity e. The values of these 
parameters for the eight models investigated in this paper 
are shown in table [lj 


3 RESULTS 

In this section, wepresent the results of the simulations for 
the fallback phasqj. The time is scaled by the period of the 
star P* with the starting point (£/P* = 0) being when the 
most bound debris of the stream come back to pericentre. 
The circularization process is investigated through the evo¬ 
lution of both the position and the specific orbital energy 
of the debris. They will be respectively compared to the so 
called circularization radius and specific circularization en¬ 
ergy. 

These two quantities are defined as follows. Before the 
disruption, the distribution of angular momentum of the gas 
inside the star is sharply peaked around that of the star. As 
this distribution is not significantly affected by the disrup¬ 
tion, the debris still have similar angular momenta. Assum¬ 
ing that they then move from their initial eccentric orbits 
to circular orbits each of them conserving their angular mo¬ 
mentum, these circular orbits will form at similar radii. It al¬ 
lows to define a characteristic circularization distance, called 
the circularization radius, obtained by equating the specific 
angular momentum of the star, given by equations J5J) and 


3 Movies of the simulations pre¬ 

sented in this paper are available at 

http://home.strw.leidenuniv.nl/~bonnerot/research.html 
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Table 1. Name and parameters of the different models. 


Model 

Potential 

EOS 

p 

e 

RI5e.8 

Relativistic 

Locally isothermal 

5 

0.8 

KI5e.8 

Keplerian 

Locally isothermal 

5 

0.8 

RA5e.8 

Relativistic 

Adiabatic 

5 

0.8 

KA5e.8 

Keplerian 

Adiabatic 

5 

0.8 

Rile.8 

Relativistic 

Locally isothermal 

1 

0.8 

RAle.8 

Relativistic 

Adiabatic 

1 

0.8 

RI5e.95 

Relativistic 

Locally isothermal 

5 

0.95 

RA5e.95 

Relativistic 

Adiabatic 

5 

0.95 



t/P* 




Figure 3. Evolution of the average specific orbital energy of the 
debris for models RI5e.8 and KI5e.8. For these models, the period 
of the star P* is 2.8 h. The average specific orbital energy is shown 
relative to the specific circularization energy given by equations 
and CE3 for the Keplerian and relativistic potentials respec¬ 
tively. 


Figure 4. Specific orbital energy distributions of the debris at 
different times t/P* = 0, 2, 4, 6 and 8 for model RI5e.8 (upper 
panel) and KI5e.8 (lower panel). For these models, the period of 
the star P* is 2.8 h. The specific orbital energy is shown relative 
to the specific circularization energy given by equations Cl and 
(I15D for the Keplerian and relativistic potentials respectively. 


evaluated at apocentre, and the specific circular angular 
momentum given by 

lc = (GM h R) 1/2 , (10) 


3.1 Impact of relativistic precession 

As mentioned in the previous section, relativistic precession 
modifies the trajectory of the debris at pericentre. This effect 


is the strongest for the orbit with /3 = 5 and e = 0.8 where 
the precession angle reaches 89.7 degrees. This is because 
this orbit has the lowest pericentre distance P p = 9.4P g and, 
to a smaller extent, the lowest eccentricity. We investigate 
the role of relativistic precession for a star on this orbit 
by comparing models RI5e.8 and KI5e.8. In model RI5e.8, 
the relativistic potential is used and apsidal precession is 
taken into account. In model KI5e.8, the debris are forced 
to move on closed orbits by using a Keplerian potential. In 
both models, a locally isothermal EOS is adopted for the 
gas. 

We discuss model RI5e.8 first. The evolution of the de¬ 
bris can be seen in Fig. [2] (upper panel), which shows snap¬ 
shots of their fallback towards the black hole at different 
times t/P* = 0,1.2, 3 and 8. The stream remains unaffected 
until t/P* ~ 1.2, corresponding to its second pericentre pas¬ 
sage. At this time, the stream self-crosses due to apsidal 
precession leading to the formation of shocks that convert 
kinetic energy into thermal energy. However, as a locally 
isothermal EOS is used, this excess thermal energy is in¬ 
stantaneously removed from the gas. The net effect of these 
shocks is therefore to reduce kinetic energy. This results in 
a decrease of the average specific orbital energy of the de¬ 
bris, whose evolution is shown in Fig. 0 (solid black line). 
As self-crossings occur at each pericentre passage, this de¬ 
crease continues for t/P* > 1.2 through periodic phases. 
These phases of decrease also become progressively steeper 
as the self-crossings involve a larger fraction of the stream. 
Accordingly, the specific orbital energy distribution, shown 
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where v a is obtained from equations J7J) or © depending on 
the potential. The associated specific orbital energy, called 
specific circularization energy, is equal to the specific circular 
orbital energy evaluated at the circularization radius. In the 
two potentials, it is given by 
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Figure 5. Snapshots of the fallback of the debris at different times t/P * = 0, 1.2, 3 and 8 for model RA5e.8. For this model, the period 
of the star P* is 2.8 h. The colours correspond to the column density £ of the gas whose value is indicated on the colour bar. The white 
point represents the black hole. The dashed white circle on the last snapshot represents the circularization radius given by equation m 
for the relativistic potential. The dotted white circle represents the semi-major axis of the star. 






^circ) / ^circ) 



Figure 6. Specific orbital energy distribution (upper panel) and 
specific angular momentum distribution (lower panel) of the de¬ 
bris at different times t/P* = 0, 2, 4, 6 and 8 for model RA5e.8. 
For this model, the period of the star P* is 2.8 h. The specific or¬ 
bital energy is shown relative to the specific circularization energy 
given by equation CESJ for the relativistic potential. The specific 
angular momentum is shown relative to that of the star given by 
equation cu evaluated at apocentre. 

in Fig. [H (upper panel), shifts towards lower energies. Ow¬ 
ing to this orbital energy decrease, the debris progressively 
move from their initial eccentric orbits to circular orbits. 
At t/P* ~ 8, their average specific orbital energy reaches a 



t/P* 

Figure 7. Evolution of the average specific orbital energy of all 
the debris, the non-accreted and accreted ones for model RA5e.8. 
For this model, the period of the star P* is 2.8 h. The average spe¬ 
cific orbital energy is shown relative to the specific circularization 
energy given by equation dj for the relativistic potential. 

value similar to the specific circularization energy (Fig. [3j 
solid black line). By this time, they have settled into a thin 
and narrow ring of radius comparable to the circularization 
radius (Fig. [2] upper panel). 

As can be seen from Figs. [2] (upper panel) and [3] (solid 
black line), the final specific orbital energy of the debris 
is somewhat larger than the specific circularization energy, 
which results in a ring forming slightly outside the circular¬ 
ization radius. These small discrepancies are due to redis¬ 
tribution of angular momentum between the debris during 
the shocks where a fraction of them (3 % at t/P* = 8) loses 
enough angular momentum to be accreted onto the black 
hole. This causes an excess angular momentum shared by 
the remaining debris, that therefore settle into circular or¬ 
bits at radii larger than the circularization radius. Notably, 
these discrepancies are reduced when the resolution of the 
simulations is increased, as will be shown in subsection 13.41 

We now discuss model KI5e.8 for which the evolution of 
the debris, shown in Fig. [2] (lower panel), is very different. As 
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Figure 8. Cross-section in the R — z plane of the torus formed at 
t/P* =8 for model RA5e.8. z represents the height with respect 
to the midplane and R the cylindrical radius. The distances are 
normalized by the tidal radius. The black hole is at the origin. 
The colours correspond to the density p of the gas whose value is 
indicated on the colour bar. 


there is no apsidal precession, the stream does not self-cross. 
Therefore, no kinetic energy is removed from the debris and 
their specific orbital energy remains constant (Fig. El dashed 
red line). Their specific orbital energy distribution also stays 
peaked around its initial value although it spreads somewhat 
due to tidal effects (Fig. [U lower panel). Consequently, the 
debris do not move to circular orbits and settle instead at 
t/P* cp 8 into an elliptical strip centred around the initial 
orbit of the star (Fig. [21 lower panel). 

The comparison of models RI5e.8 and KI5e.8 shows the 
fundamental role of apsidal precession in the circularization 
process. Our more acc urate treatment of aps idal precession 
confirms the results of lHavasaki et all ( 2013 ). 


3.2 Influence of the cooling efficiency 

The models discussed in the previous subsection (RI5e.8 and 
KI5e.8) use a locally isothermal EOS for the gas. We now 
discuss models RA5e.8 and KA5e.8, that use an adiabatic 
EOS for the gas instead. The potential used is relativistic 
for model RA5e.8 and Keplerian for model KA5e.8. The star 
is on the same orbit as before with — 5 and e = 0.8. 

We discuss model RA5e.8 first. The evolution of the de¬ 
bris as they fall back towards the black hole is shown in Fig. 
Oat different times t/P* = 0,1.2,3 and 8. The stream self¬ 
crosses at t/P* cp 1.2 and experiences shocks which convert 
kinetic energy into thermal energy. This early evolution is 
similar to model R15e.8. However, the subsequent behaviour 
differs. As the EOS is adiabatic, the thermal energy pro¬ 
duced by the shocks is not removed but kept in the debris. 
As more self-crossings occur, more thermal energy is injected 
into the stream which expands under the influence of ther¬ 
mal pressure. During the expansion, the ordered motion of 
the stream is suppressed and the debris depart from their 
initial eccentric orbits. This is achieved via a redistribution 
of their orbital parameters which can be seen in Fig. [6] by 
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Figure 9. Azimuthally-averaged ratios of the pressure (upper 
panel) and centrifugal (lower panel) forces to the gravitational 
force projected in the spherical radial direction as a function of the 
cylindrical radius R at t/P* = 8 for model RA5e.8. These ratios 
are shown for different heights z with respect to the midplane: 
\z/R t \ < 0.5, 0.5 < \z/R t \ < 1, 1 < \z/R t \ < 1.5 and 1.5 < 
\z/Rt\ < 2. The distances are normalized by the tidal radius. 


a spread in the distributions of specific orbital energy (up¬ 
per panel) and angular momentum (lower panel). The latter 
also presents a tail at large angular momenta that reaches, 
at t/Pir cp 8, the specific circular angular momentum at the 
semi-major axis of the star, indicated in the lower panel of 
Fig. [6] by a vertical black dotted line. Accordingly, the de¬ 
bris settle into a thick and extended torus with most of the 
gas located between the circularization radius and the semi¬ 
major axis of the star (Fig. 0). At this time, the majority 
of the debris (91%) are still bound to the black hole as in¬ 
dicated by their negative specific orbital energies (Fig. El 
upper panel). 

Remarkably, a significant fraction of the debris are bal- 
listically accreted onto the black hole during the formation of 
the torus. The fraction of accreted gas grows roughly linearly 
to reach 26 % at t/P* ~ 8. Fig. [7| shows the evolution of the 
average specific orbital energy of all the debris (solid black 
line), the non-accreted (dashed red line) and accreted (long 
dashed blue line) ones. When all the debris are considered, 
their orbital energy decreases due to the shocks where it is 
converted into thermal energy. However, it does not reach 
the circularization energy as part of the thermal energy is 
transferred back into kinetic energy during the expansion. 
The accreted debris are those which experienced the largest 
decrease of kinetic energy during the shocks as can be seen 
from their lower orbital energy. They also have the largest 
thermal energy which is therefore advected onto the black 
hole. The thermal energy of the torus is thus reduced due to 
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both its expansion and the accretion of the hottest debris. 
Instead of decreasing, the orbital energy of the non-accreted 
debris stays constant. This is due to the accretion of debris 
with low orbital energies which results in an excess orbital 
energy shared among the non-accreted debris. As the ac¬ 
creted debris also have the lowest angular momenta, their 
accretion leads to an increase of the angular momentum of 
the remaining debris which results in a slight shift of the spe¬ 
cific angular momentum distribution to larger values (Fig. 
El lower panel). 

A more precise analysis of the torus formed at t/P* = 8 
can be done by examining its internal configuration. Fig. [H] 
shows a cross-section in the R — z plane where z represents 
the height with respect to the midplane and R the cylindrical 
radius. It exhibits funnels around the rotation axis. It also 
contains a dense inner region at z/Rt < 0.5 and R/Rt — 1 
that corresponds to the semi-major axis of the star. This 
dense region is surrounded by a more diffuse one that ex¬ 
tends from close to the black hole to R/Rt — 2.5 and z/Rt — 
1.5. To assess the internal equilibrium of the toms, we plot 
in Fig. 0 azimuthally-averaged ratios of the pressure force 
(upper panel) and centrifugal force (lower panel) to the grav¬ 
itational force projected in the spherical radial direction as 
a function of R and for different intervals of z: \z/Rt \ <0.5, 
0.5 < \z/R t \ < 1, 1 < \z/R t \ < 1.5 and 1.5 < \z/R t \ < 2. 
The ratio of the pressure force to the gravitational force is 
roughly ( F p /F g ) r ~ 0.2 in the entire toms implying that 
it is not hydrostatically supported. However, the ratio of 
the centrifugal force to the gravitational force is larger. For 
| z/Rt | < 0.5, it presents a maximum of ( F c /F g ) r — 1 at 
R/Rt — 0.3 and decreases at larger radii. The same depen¬ 
dence exists for regions further from the midplane but for 
lower maxima of (F c /F g ) r — 0.5 — 0.8. Therefore, this toms 
is mostly centrifugally supported against gravity with this 
support being stronger in its inner region. As they are bound 
to the black hole, the regions that are not supported against 
gravity will stop expanding and collapse at a later time. This 
collapse is likely to cause the formation shocks in the outer 
part of the toms which would increase hydrostatic support 
in this region. 

For model KA5e.8, the stream does not experience apsi- 
dal precession. However, it is heated when it passes at peri- 
centre. This is due to the formation of the pancake shock de¬ 
scribed in the Introduction. As a result, it expands roughly 
by a factor of 2 at each pericentre passage both in the radial 
and vertical directions. It causes the stream to self-cross at 
t/P * ~ 4 that is after its fifth passage at pericentre. The 
subsequent evolution is similar to model RA5e.8. The de¬ 
bris settle into a thick and extended toms at t/P* ~ 10. 
This channel of disc formation is similar to that found by 
iRamirez-Ruiz &; RosswoeJ (l2009l h Therefore, relativistic ap- 
sidal precession is not the only factor that can lead to cir¬ 
cular orbits of the debris, but it is the most efficient and 
operates regardless of the EOS used for the gas. 

As can be seen by comparing models RI5e.8 and RA5.8, 
the cooling efficiency determines the structure of the disc 
formed during the circularization process. While a thin and 
narrow ring forms at the location of the circularization ra¬ 
dius for an efficient cooling, an inefficient cooling leads in¬ 
stead to the formation of a thick and extended toms located 
between the circularization radius and the semi-major axis 
of the star. 


3.3 Dependence on the orbit of the star 

So far, the fallback of the debris has been investigated for 
a fixed orbit of the star with f3 — 5 and e = 0.8. We now 
consider two new orbits obtained by decreasing the penetra¬ 
tion factor (subsection 13.3.II) and increasing the eccentricity 
(subsection 13.3.21) . Throughout this section, the potential 
is fixed to relativistic. Models RI5e.8 and RAe.8, discussed 
above for the initial orbit and this potential, will be used as 
reference. 


3.3.1 Decreasing the penetration factor 

We start by decreasing the penetration factor to /3 = 1 keep¬ 
ing the eccentricity to its initial value e = 0.8. This decrease 
of the penetration factor corresponds to an increase of the 
pericentre distance from R p — 9AR g to R p = 47 R g . We 
investigate the fallback of the debris for a star on this new 
orbit by discussing models Rile.8 and RAle.8. A locally 
isothermal EOS is used for model Rile.8 while an adiabatic 
EOS is used for model RAle.8. The relativistic potential is 
used in both models. 

We discuss model Rile.8 first. Snapshots of the evolu¬ 
tion of the debris are shown in Fig. [TU] at different times 
t/P* = 0, 2.2,4, 6,8 and 16. The first self-crossing of the 
stream occurs at t/P * = 2.2 that is at its third pericentre 
passage. It causes the formation of shocks that reduce the 
kinetic energy of the debris as a locally isothermal EOS is 
used. However, as the stream passes further from the black 
hole, relativistic precession is weaker than for model RI5e.8 
with a precession angle decreasing from 89.7 to 13.5 degrees. 
Consequently, the shocks happen further from the black hole 
and involve parts of the stream with lower relative velocities. 
For this reason, they are less efficient at removing kinetic 
energy from the debris. This results in a slower and more 
gradual decrease of their average specific orbital energy, as 
seen from Fig. [IT] that shows its evolution for model RI5e.8 
(solid black line) and Rile.8 (dashed red line). At t/P* ~ 16, 
the average specific orbital energy of the debris stabilizes at 
a value similar to the circularization energy (Fig. [IT] dashed 
red line) and the debris settle into a thin and narrow ring 
of radius comparable to the circularization radius (Fig. [TUT). 
As for model RI5e.8, the origin of the small discrepancies 
with the specific circularization energy and the circulariza¬ 
tion radius are due to the accretion onto the black hole of a 
fraction of the debris (3 % at t/P* — 16) with low angular 
momenta. 

For model RAle.8, the evolution of the debris is very 
similar to model RA5.8. The shocks produced by the self¬ 
crossings of the stream leads to its expansion. The debris 
settle into a torus with most of them located between the 
circularization radius and the semi-major axis of the star. 
The only effect of modifying the orbit is to change the value 
of these two limiting radii. 


3.3.2 Increasing the eccentricity 

We now keep the penetration factor to its initial value /3 = 5 
while increasing the eccentricity to e = 0.95. The fallback 
of the debris is investigated for this new orbit by discussing 
models RI5e.95 and RA5e.95 which use respectively a locally 


© 2014 RAS, MNRAS OOO.fTHTTI 








10 Clement Bonnerot, Elena M. Rossi, Giuseppe Lodato and Daniel J. Price 



5 6 7 

log 2 (g/cm 2 ) 


Figure 10. Snapshots of the fallback of the debris at different times t/P* = 0, 2.2,4, 6,8 and 16 for models Rile.8. For this model, the 
period of the star P* is 31 h. The colours correspond to the column density E of the gas whose value is indicated on the colour bar. The 
white point represents the black hole. The dashed white circle on the last snapshot represents the circularization radius given by equation 
03 for the relativistic potential. Apsidal precession is weaker for this model than for model RI5e.8 owing to a larger pericentre, which 
causes the stream to self-cross further out from the black hole. At this location, relative velocities are lower, weakening the shocks. The 
debris therefore move slower to circular orbits. 


isothermal and an adiabatic EOS. The relativistic potential 
is adopted for both of them. 

We discuss first model RI5e.95. Fig. [12] shows snap¬ 
shots of the fallback of the debris at different times t/P* = 
0,0.15,0.25,0.33,0.5 and 1. The stream first crosses itself 
at t/P* ~ 0.15 just after the first passage of its leading 
part at pericentre. When this occurs, most of the debris are 
still falling back towards the black hole away from the self¬ 
crossing point. This is different from model RI5e.8, for which 
this self-crossing happens later and when most of the debris 
are already beyond the self-crossing point (see Fig. [2]). The 
reason for these differences is that the stream is longer for 
model RI5e.95 than for model RI5e.8 owing to the larger 
eccentricity. After the first self-crossing, the debris located 
beyond the self-crossing point are expelled leaving the rest of 
the stream free to move around the black hole. This induces 
a second self-crossing at t/P* ~ 0.33. These self-crossings 
create shocks that remove kinetic energy from the debris as 



Figure 11. Evolution of the average specific orbital energy of the 
debris for models RI5e.8, Rile.8 and RI5e.95. For these models, 
the periods of the star P* are respectively 2.8 h, 31 h and 22 h. 
The average specific orbital energy is shown relative to the specific 
circularization energy given by equation m for the relativistic 
potential. 
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Figure 12. Snapshots of the fallback of the debris at different times t/P* = 0,0.15,0.5 and 1 for models RI5e.95. For this model, the 
period of the star P* is 22 h. The colours correspond to the column density S of the gas whose value is indicated on the colour bar. 
The white point represents the black hole. The dashed white circle on the last snapshot represents the circularization radius given by 
equation for the relativistic potential. The stream is longer for this model than for model RI5e.8 as the eccentricity is larger. As a 
result, the self-crossings involve a larger fraction of the debris, causing them to move faster to circular orbits. 


a locally isothermal EOS is used. As shown in Fig. [IT] (blue 
long dashed line), the average specific orbital energy of the 
debris therefore decreases for t/P* > 0.15. This decrease oc¬ 
curs faster than for model RI5e.8 as the self-crossings involve 
larger fractions of the debris. At t/P* ~ 1, their average spe¬ 
cific orbital energy settles at a value similar to the circular¬ 
ization energy (Fig.[TTJ blue long dashed line) and they form 
a thin and narrow ring of gas located at a distance compara¬ 
ble to the circularization radius iFig.1121 ) . As previously, the 
discrepancies with the specific circularization energy and the 
circularization radius are due to the accretion of a fraction 


of the debris (4 % at t/P* = 1) with low angular momenta. 
These discrepancies are larger for this model than for mod¬ 
els RI5e.8 and Rile.8 because more debris are accreted and 
these debris have lower angular momenta relative to that of 
the star. 

For model RA5e.95, the debris evolve very similarly to 
model RA5.8. They expand due to shocks produced by the 
self-crossings of the stream. They settle into a thick and 
extended torus in which most of the debris are located be¬ 
tween the circularization radius and the semi-major axis of 
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the star. Only the value of these two limiting radii differs 
for the new orbit. 

3.4 Convergence of the results 

The simulations of the fallback phase has been performed 
for three different resolutions corresponding to about 100K, 
500K and 1300K particles. The two larger resolutions have 
been obtained fr om the lowest on e by using the particle split¬ 
ting technique (|Kitsionasl 12003 ) on the initial condition of 
the fallback phase. Each particle is split into 5 or 13 addi¬ 
tional ones with one at the position of the initial particle 
and the others distributed on a tetrahedron and on a face- 
centred cubic structure respectively. In both cases, the dis¬ 
tance between each additional particles is fixed to 1.5 times 
the smoothing length of the initial one. In addition, the par¬ 
ticles added outside the volume defined by the initial ones 
have been removed. The fraction of particles removed in this 
way is always less than 10%. The mass of the additional par¬ 
ticles has then been decreased in order to keep the total mass 
of the stream constant. 

Minor differences have been noticed when the resolu¬ 
tion is increased. They are common to the different orbits 
considered. For a locally isothermal EOS, the time at which 
the average specific orbital energy of the debris settles is 
either advanced or delayed by At ~ P* between the 100K 
and the 500K simulations. However, this difference becomes 
negligible between the 500K and 1300K simulations. Both 
for a locally isothermal and an adiabatic EOS, increasing 
the resolution also results in a smaller fraction of debris 
with low angular momenta being accreted onto the black 
hole. For the isothermal EOS, this causes the average 
specific orbital energy of the debris to settle at a value 
closer to the specific circularization energy. The fraction 
of debris accreted during the disc formation decreases to 
less than 5% in this case. We therefore attribute it to 
resolution. For the adiabatic EOS, the torus turns out to be 
slightly more compact and centrally condensed. However, 
its internal structure remains the same. The fraction of 
accreted debris also decreases in this case but remains larger 
than 25%, and is likely to be physical. The simulations 
presented in this paper have been performed for 500K 
particles as both the final configuration of the debris and 
the time needed to reach it are unchanged above this reso¬ 
lution. Overall, we conclude that the behaviour described in 
the previous subsections is robust with respect to resolution. 


4 DISCUSSION AND CONCLUSION 

By means of SPH simulations, we investigated the circular¬ 
ization process for tidal disruptions of stars on bound orbits 
by a non-rotating black hole. The formation of an accretion 
disc from the debris is mostly driven by relativistic apsidal 
precession that causes the stream to self-cross. If cooling 
is inefficient, this self-crossing is also partially caused by an 
expansion of the stream, when it passes at pericentre. In ad¬ 
dition, we showed that the structure of the disc depends on 
the cooling efficiency of the gas by considering two extreme 
cases. For an efficient cooling, the debris form a thin and nar¬ 
row ring of gas. For an inefficient cooling, they settle into a 


thick and extended torus, that at the end of our simulation is 
still mostly centrifugally supported against gravity. We also 
demonstrated the existence of different regimes of circular¬ 
ization for different orbits of the star. The circularization 
timescale t c ire varies from ~ P* for the largest eccentricity 
to ~ 10P* for the lowest penetration factor considered (see 
Fig. Hill. In physical units, it corresponds to tens to hundreds 
of hours. 

The circularization proces s in T DEs has been th e focus 
of several other recent works. IShiokawa et al.l (120151 ) simu¬ 
lated the tidal disruption of a white dwarf on a parabolic 
orbit around a 500 M 0 black hole us i ng a g eneral relativis¬ 
tic simulation while lHavasaki et al.l (120151 ) considered the 
case of bound orbits around a 10 6 M 0 black hole. The lat¬ 
ter use a simulation setup similar to ours but our treatment 
of apsidal precession is more accurate and we consider a 
larger range of orbits. Both studies noticed the influence 
of the cooling efficien cy on the disc structure. In addition, 
lHavasaki et al.1 (|2015l ) invest igated the effect of the black 
hole spin. As anticipated by I Dai et af ] d2013h . they found 
that nodal precession can prevent the self-crossing of the 
stream, delaying circularization. Us i ng Mo nte-Carlo calcu¬ 
lations, iGuillochon & Ramirez-Ruia (120151 ) found that this 
delay is of one year on average for a star on a parabolic orbit. 

4.1 Thermal energy radiation 

During the circularization process, thermal energy is in¬ 
jected into the debris at the expense of their kinetic en¬ 
ergy. If they cool efficiently, this thermal energy increase 
can be estimated by subtracting the final circularization en¬ 
ergy to the initial orbital energy of the star. It amounts to 
A U = — £^ rc ) — few 10 51 — 10 52 erg for the orbits 

considered, according to equations and csi). For larger 
A U is larger because the circularization energy is lower. 
Assuming that this excess thermal energy is radiated dur¬ 
ing the disc formation, it gives rise to a flare of luminosity 
L C irc = AU/ t c i rc ~ few 10 — 10 3 LecM, increasing with /3, as 
tcirc is shorter and AU is larger. 

As described in subsection cm if the excess thermal 
energy is not radiated instantaneously, it is partially trans¬ 
ferred back into kinetic energy or advected onto the black 
hole. As a result, the remaining thermal energy of the torus 
is ~ 10% of the above AU. The associated heating rate is 
therefore AU/Girc — few 1 — 10 2 LEdd- Note that at the end 
of our simulation, the torus has not yet settled into an equi¬ 
librium configuration as we observe outflowing material. 

4.2 Viscous evolution 

Debris on circular orbits are subject to viscous effects which 
drive their accretion onto the black hole. Viscosity is not 
explicitly included in the simulations presented above. Nev¬ 
ertheless, its influence can be estimated a posteriori by com¬ 
puting the viscous timescale Giso In our simulations, disc 
formation is progressive. A fraction of the debris already 
forms a disc-like structure before complete circularization. 
An example of such a structure can be seen in Fig. [TJ] at 
£/P* = 0.33. We consider that viscosity can start acting on 
these debris at t < t c ire which may lead to their accretion 
before complete circularization if t v isc < Giro For simplic¬ 
ity, tvisc is computed for the final configuration of the disc 
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although it is used to evaluate viscous effects as the debris 
proceed towards this configuration. If the debris cool effi¬ 
ciently, they settle into a thin ring around the circularization 
radius. At this distance from the black hole, 


tt= 3xio, (< &r (pH 


1-e 2 

0.36 


3/2 


(16) 


wher e a is the viscosity parameter ([Shakura &; Sunvaevl 
1973) and H/R is the aspect ratio of the disc. The circu¬ 
larization radius is obtained from equation © which is at 
the origin of the (1 — e 2 ) 3//2 factor. For all orbits consid¬ 
ered, tvisc > tcirc, largely independently of the values of a 
and H/R. Therefore, viscosity should not have a significant 
effect during circularization and most of the accretion will 
occur once the ring is formed. At the end of circularization, 
the viscous timescale, given by equation G9, corresponds 
to an accretion rate M = M*/t V isc — few 1 — lOOMEdd de¬ 
pending on the orbit, MEdd being the Eddington accretion 
rate for a radiative efficiency of 10%. We therefore specu¬ 
late that, due to this super-Eddington accretion, the ring 
will subsequently evolve into a thicker structure under the 
influence of radiation pressure. 

If the gas cools inefficiently, it settles into a thick torus 
located between the circularization radius and the semi¬ 
major axis of the star. At the semi-major axis of the star, 




(17) 


which does not depend explicitly on the orbit of the star 
as P* cancels out. For the above values of a and H/R , 
t v isc < Pirc which is sensitive to the values chosen for these 
parameters. However, when evaluated at the circularization 
radius, t v isc is lower by a factor (1 — e 2 ) 3 / 2 ~ few 10 2 —10 1 
depending on the orbit. The precise value of the viscous 
timescale thus depends on the mass distribution within the 
torus, being shorter if more mass is close to the circulariza¬ 
tion radius. We therefore conclude that viscosity may affect 
the evolution of the debris during the circularization process 
causing some of them to be accreted. During its subsequent 
evolution, the torus will keep accreting matter with an ac¬ 
cretion rate M = M*/£ v isc > 10 4 M^dd for the different 
orbits considered. For these highly super-Eddington accre¬ 
tion rates, the subsequent evolution of the torus is difficult 
to predict. 


4.3 Evaluation of the cooling efficiency 

In the simulations presented in this paper, two extreme cool¬ 
ing efficiencies have been considered. The ability of the de¬ 
bris to cool can be estimated a posteriori by computing the 
diffusion timescale tdiff defined as the time that photons take 
to diffuse out of the surrounding gas. As most of the ther¬ 
mal energy is produced by the shocks occurring when the 
stream self-crosses, the diffusion timescale must be evalu¬ 
ated at the location of these shocks. Since we find that the 
gas is optically thick to electron scattering, it is given by 
^diff = H sh r/c where r = cr T p sh P s h/m p is the optical depth, 
H s h and p s h are the width and density of the gas at the lo¬ 
cation of the shocks. When an efficient cooling is assumed, 
the width of the stream remains H s h — R*. To ensure self- 


consistency, the condition tdis < £ c irc has to be satisfied, 
which translates into an upper limit on the density 


Psh <8x10 


cm 




a* \ 

100 R©y 


3/2 


(18) 


where n C irc = tcirc/P*- Among the different orbits consid¬ 
ered, this condition is satisfied only for that with [3 — 1. 
In general, we estimate that increasing the eccentricity e 
of the star tends to favour an efficient cooling of the de¬ 
bris. This is because it leads to a more extended and tenu¬ 
ous stream, decreasing p s h. Condition (TTB1) is therefore more 
easily fulfilled. Instead, increasing the penetration factor [3 
may favour an inefficient cooling. It causes the stream to 
self-cross closer to the black hole. The debris are therefore 
located in a smaller volume which increases p s h. Further¬ 
more, they circularize faster which decreases n C i rc - Condition 
m is thus more difficult to fulfill. However, these estimates 
are only approximate. The ability of radiation to escape de¬ 
pends on the precise location within the stream at which 
the thermal energy is deposited by the shocks. It also varies 
with the density distribution of the debris which is highly 
inhomogeneous at the time when most of the shocks occur. 
Furthermore, this radiation could also affect the structure of 
the disc through radiation pressure. A realistic treatment of 
the interaction between gas and radiation is therefore nec¬ 
essary to determine precisely the influence of cooling during 
the circularization process. 


4.4 Extrapolation to parabolic orbits 

As mentioned in the Introduction, TDEs typically involve 
stars on parabolic orbits. In this paper, we chose to simulate 
the disruption of stars on bound orbits instead. This choice 
is motivated by a lower computational cost which allows 
to explore a larger parameter space. However, the results 
can be extrapolated to get insight into the typical case of 
parabolic orbits. 

The bound orbits considered in the simulations satisfy 
the condition e < e cr it = 1 — (2//3)(Mh/M*) _1//3 , so that all 
the debris are bound to the black hole. Furthermore, all the 
debris have periods similar to that of the star, P*. Instead, 
in the case of parabolic orbits, part of the debris become 
unbound from the black hole. Therefore, the tidal stream 
has a large range of periods between that of the most bound 
debris t m i n and +oo. This means that it features an elon¬ 
gated tail of debris which will continue to fall back towards 
pericentre long after the most bound ones have reached it. 
Due to apsidal precession, the leading part of the stream 
will inevitably collide with this tail after its first passage 
at pericentre. By means of point particle calculations in 
the Schwarzschild metric, we found that this prompt self- 
crossing occurs in general for an eccentricity e > 0.9 largely 
independent on (3. As discussed in subsection 13.3.21 it hap¬ 
pens for model RI5e.95 where e = 0.95 and /3 = 5 (see Fig. 
rm In this case, the self-crossing leads to circularization on 
a timescale t C irc — P*. Therefore, for a parabolic orbit with 
/3 = 5, we expect a disc to form from the most bound debris 
due to this prompt self-crossing on a timescale t c ire — tmin 
where t m i n replaces P* as the period of the most bound de¬ 
bris. In this case, the newly-formed disc is only composed 
of the tip of the stream, which falls back within ~ t m in- We 
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expect the debris infalling later to rapidly circularize and 
join this disc. The timescale t c ire for the most bound debris 
to circularize is less clear for a parabolic orbit with /3 = 1 as 
apsidal precession is weaker. For model Rile.8 where e = 0.8 
and /3 = 1, it is t C irc — 10P*. However, the presence of a tail 
of debris would cause the self-crossing to affect the middle of 
the stream instead of its extremities, which likely makes the 
shocks more disruptive. On the other hand, the self-crossing 
happens further out from the black hole. This is because the 
apocentre of the most bound debris is about ten times larger 
for a parabolic orbit than for a bound one with e = 0 . 8 . 
Relative velocities are lower at this location, which likely 
weakens the shocks. Therefore, t c ire may remain ~ 10£ m in 
for a parabolic orbit with /3 = 1. In this case, a significant 
mass of debris falls back within t c ire, while the most bound 
ones circularize. Whether all this mass has circularized by 
tdre is unclear. 

Based on our estimates in subsection 14.31 conditions for 
efficient cooling are more easily met for parabolic orbits than 
for elliptical ones. However, how efficient the cooling is also 
depends on the penetration factor, as large /3 favours the 
formation of denser regions. In general, we expect the cooling 
efficiency to have the same effect on the disc structure as 
found in our simulations. If cooling is efficient throughout 
the evolution, a thin ring forms around the circularization 
radius. If it is inefficient, a thick torus forms located between 
the circularization radius and around the semi-major axis of 
the most bound debris. This allows us to extrapolate to the 
case of parabolic orbits the discussion in subsections 14.11 and 

roi 

If cooling is efficient, the increase of thermal energy ex¬ 
perienced by the debris during circularization is equal to 
the total change in orbital energy of the most bound de¬ 
bris. It amounts to A U — /M*(£ R — £ R rc ) — few 10 51 erg 
according to equations ® and (USD For (3 = 1, this esti¬ 
mate assumes that all the debris falling back within t C irc 
have circularized. The factor / accounts for the fact that 
only a fraction of the debris reaches the black hole within 
tcirc in the parabolic case. Numerically, / ~ 0.2 — 0.4 for 
tdre = 1 — 10 tmin assuming a flat energy distribution. If this 
thermal energy is radiated during t c ire, it leads to a lumi¬ 
nosity Z/ C irc = AP/tcirc — few 1 — lOLEdd- Remarkably, it is 
comparable or higher than the peak luminosity in the soft 
X-ra y band from the viscou s accretion accretion through the 
disc (iLodato &; Rossilfeoilh . If the thermal energy is not im¬ 
mediately radiated but partly used to expand the disc or 
advected onto the black hole, the remaining thermal energy 
of the toms is ~ 10 % of the above A U, a percentage ex¬ 
trapolated from the case of bound orbits. The heating rate 
is then AU/t C i rc — fe w 0.1 — lLsdd- This ro ughly agrees 
with the value found bv IShiokawa et all (l2015h scaling their 
results to ou r parameters for the black hole and the star 
(jPiran et al.ll2015h . 

The effect of viscosity can be estimated in the case of 
parabolic orbits by computing the viscous timescale p isc, 
obtained from equations (fT(Jl) and (1171) replacing P* and e 
by the period and eccentricity of the most bound debris t m in 
and e cr it- As for bound orbits, we expect the disc to form 
progressively, with a fraction of the debris rapidly reach¬ 
ing a disc-like structure. For ft = 1, these circularized de¬ 
bris could already be present before most of the final disc 
mass reaches pericentre. On these debris, viscosity may start 


acting even before the disc completely settles. If cooling is 
inefficient, we find that pisc — few t m i n at the semi-major 
axis of the most bound debris but decreases by a factor 
(1 — e^it ) 3 / 2 — few 10 -3 — 10 -4 at the circularization radius. 
The distribution of mass within the torus thus determines 
the value of the relevant viscous timescale which decreases 
if most of the debris are close to the circularization radius. 
In all cases, we expect pisc < Pirc since t c ire is estimated 
to be > tmin • We conclude that a significant fraction of the 
toms may already be accreted by the end of circularization. 
If cooling is efficient, pisc — few 10 — 100 t m i n in the thin 
ring. In this case, pi SC > Pire which indicates that the ring 
is not significantly accreted during circularization. However, 
if the disc geometry changes when accretion starts, the vis¬ 
cous timescale could be different. In particular a thicker disc 
may result from a super-Eddington accretion rate, shorten¬ 
ing the viscous timescale. 

The fact that pi SC tmin has often been interpreted as 
evidence that the accretion rate onto the black hole traces 
directly the fallback rate of the debris. However, if most of 
the debris initially form a disc without having the possibility 
of accreting, then the accretion rate onto the black hole will 
be driven by viscous processes rather than by infall, which 
has essentially finished once accretion starts (| Cannizzo et al.l 
ll990l : IShen fc Matznerll2014h . If disc accretion is significantly 
delayed wit h respect to disc form ation, we thus expect a 
solution a la lCannizzo et all (|l990h . that predicts an accre¬ 
tion rate declining as If on the other hand accretion 

already occurs while the disc is forming, we expect at late 
times a t~ 5 ^ 3 decline in the accretion rate. Owing to a pro¬ 
gressive disc formation, we consider the latter scenario to be 
a possibility, which seems favoured if cooling is inefficient, 

Since tvisc P ire* 

The simulations presented in this paper allowed to get 
insight into the circularization process during TDEs. Sev¬ 
eral extensions of this work are possible including simula¬ 
tions considering a more eccentric star, the treatment of the 
interaction between the debris and the radiation emitted 
during the disc formation and the effect of the black hole 
spin. 
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